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Abstract. We study transport in ferromagnetic single-electron transistors. The non-equilibrium spin ac- 
cumulation on the island caused by a finite current through the system is described by a generalized 
theory of the Coulomb blockade. It enhances the tunnel magnetoresistance and has a drastic effect on the 
time-dependent transport properties. A transient decay of the spin accumulation may reverse the electric 
current on time scales of the order of the spin-flip relaxation time. This can be used as an experimental 
signature of the non-equilibrium spin accumulation. 

PACS. 73.40. Gk Tunneling - 75.70.-i Magnetic films and multilayers - 73.23.Hk Coulomb blockade; 
single-electron tunneling - 75. 70. Pa Giant magnetoresistance 



■ 1 Introduction the capacitance of the island), and the tunnel resistances 



Single-electron tunneling has been an active area of re- 
search during the last decade (for a review see Ref. p|). 
In a double tunnel junction system, the Coulomb blockade 
effect is pronounced when the charging energy of the is- 
land is larger than the thermal energy, e 2 /2C > fcsT (C is 
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are larger than the quantum resistance R > Rk = h/e 2 . 
So far most of the research on the double tunnel systems 
has been in systems using metals (normal or supercon- 
ducting) or semiconductor hetero-structures Jt],^]. The so- 
called orthodox theory of single-electron tunneling jjj has 
been very successful in explaining experiments where the 
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single-particle energy separation is smaller than the ther- have also been seen in discontinuous multilayers || and in 

mal energy. small cobalt clusters [fl0|| . 

Quite unrelated, the giant magnetoresistance in mag- A generalization of the orthodox theory to describe a 

netic multilayers and the spin-tunneling magnetoresistance FSET by introducing spin-dependent tunneling rates has 
in ferromagnet-insulator-ferromagnet systems have attractedbeen presented by Barnas and Fert [|ll| and Majumdar 

great interest ■ In a magnetic multilayer or in a tri- and Hershfield [Q . However they neglected the effects of 

layer "spin- valve" structure the magnetization of the ad- a non-equilibrium spin accumulation on the island caused 

jacent magnetic layers may vary. When the layers are an- by the spin-dependent tunneling rates, which can have a 

tiparallel, an external magnetic field can align the magne- drastic effect on the transport properties of the ferromag- 

tization of the ferromagnetic layers. The observed decrease netic single-electron transistor. 

in the resistance is a result of spin-dependent scattering T r ,11 , , 

^ ^ b In a ferromagnetic smgle-electron transistor, the tran- 

in the materials. . . c t c ^ ■ 1 j j 

sition rates tor tunneling into or out of the island de- 

In recent experiments Ono et al. studied the proper- pend on the electron spin. An electric current through 
ties of ferromagnetic single-electron transistors (FSETs) the system therefore implies a spin current out of or into 
|,|. This is a new and interesting system in which both the island (ds/dt)^. This creates a non-equilibrium ex- 
the charging aspect due to the Coulomb energy of the cess spin on the island. The excess spin s decays by spin- 
excess charge on the island and the magnetoresistance flip relaxation so that s = (ds/dt)^ r T s ^, where r g r is the 
due to the spin-dependent tunneling rates are of impor- spin-flip relaxation time. The non-equilibrium spin accu- 
tance. Ono et al. found an enhanced magnetoresistance mulation on the island is equivalent to a non-equilibrium 
in the Coulomb Blockade regime and a monotonic phase chemical potential difference Afi between the spin-up and 
shift of the Coulomb oscillations induced by the mag- the spin-down states. In the case of a normal metal is- 
netic field. The enhancement of the magnetic valve ef- land this chemical potential difference can be evaluated 
feet in the Coulomb blockade regime has been ascribed as A/i = sS, where S is the single-particle energy spac- 
to co-tunneling, which is of a higher order in the tunnel- ing at the Fermi energy (inverse density of states). The 
ing resistances and so the difference in the spin-dependent chemical potential difference modifies the transport prop- 
tunneling between the antiparallel and the parallel config- erties. The non-equilibrium spin accumulation is impor- 
uration is larger ||[7||. The magneto-Coulomb oscillations tant for the transport properties of the FSET when the 
can be explained in terms of changes in the free energy non-equilibrium chemical potential difference A\i is of the 
of the island electrode and the leads in the presence of same order as the Coulomb charging energy E c = e 2 /2C. 
an external magnetic field M . Coulomb charging effects The spin current is of the same order as the electric cur- 
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rent, e(ds/dt)^ r ~ /. We are interested in voltages of the for a large island, our results reduce to the models in Refs. 
order of the Coulomb gap, so we have I ~ GE c /e, where j[l],[l^], where the spin was assumed to be equilibrated. 
G is the typical junction conductance. From these crude Spin-accumulation is important in small systems also 
considerations we deduce that the non-equilibrium spin when the boundary roughness gives a dominant contribu- 
accumulation is important when the spin-flip relaxation tion to the spin-flip scattering. In such a case, the spin- 
time and/or the single-particle energy separation are suf- orbit relaxation time is rso = /t, where / is the ratio of 
ficiently large, i.e. when non spin-flip to spin-flip scattering rates and r is the mo- 
mentum relaxation time. For a small system with an ide- 



ally rough surface r ~ L/vp, where L is the size of the sys- 
tem and vf is the Fermi velocity. Abrikosov and Gor'kov 
calculated g§ / = {aZ)~ 4 , where a = e 2 /(hc) ~ 1/137 is 
the fine-strucure constant, in rough agreement with exper- 



iments 1 17 . Using L ~ aiV 1 / 3 , where a ~ 1/kp is the lat- 
tice constant r ~ N 1 / 3 / {TiEp) and the relation ([!]) reads 
for small particles with rough boundaries: 



r s{ S/h > R/R K , (1) 

irrespective of E c . 

The spin-flip relaxation time is crucial for the observa- 
tion of the non-equilibrium spin accumulation. The spin- 
flip scattering in clean metals at low temperatures is dom- 
inated by the spin-orbit scattering, r s f w Tgo- Since spin- 
orbit coupling is a relativisitic effect, it increases with the 

atomic number Z of the metal. The single-particle energy ^ ^ R J ' ^ 
spacing on the island increases with decreasing number Assuming that the tunnel resistance is of the order of the 
of atoms N on the island as 5 ~ Ep/N, where Ep is the quantum resistance (e.g. in the measurement by Ono et 
Fermi energy. Thus the size of the island must be small for al. || R/Rr ~ l)j and using estimates for the ratio of 
heavy elements in order to observe big effects of the spin- non spin-flip to spin-flip scattering rates from the mea- 
accumulation. Let us first estimate the spin- flip relaxation surements in Ref. [17j], we expect to see effects of a non- 
time in an island with an ideal clean surface, where the equilibrium spin accumulation in Cu (log/ ~ 2.1 — 2.4) 
spin-flip scattering at the boundary of the island can be when N < 10 3 , in Al (log/ - 2.8 - 3.6) when N < 
disregarded. The spin-flip relaxation time in single-crystal 10 4 — 10 5 , in Na (log / ~ 4.2 — 5.3) when N < 10 6 — 10 9 
aluminium is (r gf - lCT 8 s at T = 4.3K fU ( T gf ~ l(r 10 s and in Li (log/ - 4.9 - 7) when A^ < 10 7 — 10 10 . For 
in polycrystalline aluminium and r g f ~ 10 _11 s for smooth boundaries, these numbers should be correspond- 
gold|l^]. The Fermi energy is Ep ~ lOeV, and the spin ingly larger. For our purposes, the study of effects of a 
accumulation may therefore be expected to play a signif- non-equilibrium spin accumulation, one should choose a 
icant role in an Al island with less than 10 s atoms (10 6 light metal with a long spin-orbit relaxation time (a large 
atoms in polycrsalline aluminium). In the limit r g r — > or parameter /). 
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Note that the orthodox model is still valid for the single-electron transistor by generalizing the orthodox the- 

small clusters described above. For a small island, the ca- ory. The simple case of a halfmetallic ferromagnet-normal 

pacitance scales like C ~ aN 1 ^ 3 . The ratio between the metal-halfmetallic Ferromagnetic (HF/N/HF) FSET can 

single-particle energy spacing and the Coulomb charging be treated analytically and is considered in section |j. The 

energy is therefore S/Ec ~ (Epa/e 2 ) N~ 2 / 3 . The prefac- transport properties in the general case are investigated 

tor is [Epa/e 2 ^ ~ 1, and thus the single-particle spacing is numerically in section ^. Finally our conclusions are pre- 

much smaller than the Coulomb charging energy. Our gen- sented in section ||. Selected results of the present work 

eralized orthodox model is therefore valid when the tern- have been presented in Ref. pp| . 
perature is larger than the single-particle energy spacing. 

For lower temperatures, the single-particle energy spacing 2 Model 
will also appear in the current-voltage characteristics. 

We consider the situation shown in Fig. |l| A (ferromag- 

" Modern" metals, like arm-chair nanotubes |L8fl or (mag- . . 

netic or normal) metal island is attached to two (ferro- 

netic) semiconductor heterostructures JL9| can also be in- 

1 — ' magnetic or normal) leads by two tunnel junctions. There 

teresting as island materials. The first because of a possi- 

is an applied source-drain voltage V between the right and 

ble huge spin-flip relaxation time (Z — 6 for carbon) and 

the left reservoir and a gate voltage source coupled capac- 

the latter since islands containing a small number of elec- 

itively to the island. We assume collinear magnetization 

trons can be created by depletion of the two-dimensional 

in the leads and the island. The direction of the magne- 

electron gas [El. However, in these systems quantum size . . . . 

LJ tizations m the island and the right lead can be parallel 

effects start to play a role which have our attention. . . .... 

or antiparallel to the direction of the magnetization m the 

In small tunneling systems where Eq. (0) is satisfied, left reservoir as depicted by the broken arrows in the Fig. 

the spin-flip relaxation time is longer than the charge re- when the coercive fields of the magne ts are different 

laxation time RC. This can be seen from Eq. @, r gf > diffcrent configurations of the magnetizations can be re- 

(2E C /8)RC, and noting that the charging energy is larger alized by sweeping an external magnetic field. We ignore 

than the single-particle energy spacing for all but in few- for simplicity the complications due to the possible ap- 

clcctron systems. Hence the long-time response of the sys- pearanC e of a magnetic domain structure in the ferromag- 

tem is dominated by the spin dynamics. We will demon- nets _ Tne tunnel j unct ions have a capacitance C, and spin- 

strate that this feature can be used as an experimental dependent conductances G ia , where i = 1,2 denotes the 

signature of the non-equilibrium spin accumulation. first and the gecond j unct ion and a denotes up (+) or down 

The paper is organized in the following way. In the (—) spin electrons. The spin-dependent conductance of the 

next section H we model the transport in the ferromagnetic first junction, Gi a , is proportional to the spin-dependent 
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density of states in the left reservoir at the Fermi level faster than the time between the tunneling events ~ RC, 
Nt (0), the island at the Fermi level iVj (0) and the angu- so that the distributions of the energy levels for the dif- 
lar averaged spin-dependent tunneling probability at the ferent spins are described by Fermi functions. The non- 
Fermi level, G\ a — 2ite 2 NL a -(0)Ni lT (0)Ti a (0)/h In a equilibrium chemical potentials in the island are spin-split 
ferromagnet the tunneling probability is significantly dif- by Afi due to the spin accumulation. The non-equilibrium 
ferent for the s- and (i-electrons. The density of states and chemical potential difference causes novel phenomena on 
the tunneling probabilities should therefore be regarded as the single-electron transistor to be discussed below, 
phenomenological parameters which are to be constant in The tunnel conductances are taken to be much smaller 
an interval of the order of the applied source-drain voltage than the quantum conductance G i<T -C Gk — I/Rk- Thus 
at the Fermi surface. This assumption has been found to we disregard co-tunneling JjJ and calculate the current- 
be valid if the applied source-drain voltage is lower than voltage characteristics in the system by using the semiclas- 
100 mV [pT[ . The expression for the spin-dependent con- sical master equation. The transition rates can be found 
ductance of the second junction is similar. We take the from Fermi's Golden Rule [[!). The spin-dependent tunnel- 
gate capacitance C g to be small compared to the junction ing rates are 
capacitances C\ and C 2 and other impedances are disre- 
garded. The energy difference associated with the tunnel- 



r w w +i,n = \G la F(E 1 {V,q)-aAn/2), (3) 



1 



ing of one electron into the island through junction i is "'" T1 e 2 

1 



r iCT „,„+i = ^G lry F(-E 1 (V,q) + crA f i/2), 



r^ n+ i, n = ^G 2 „F{E 2 (-V,q) - *An/2), 



e 



El (V,q) = Kl eV+ e{q ~ e/2 \ ^V"+i = ±G 2a F(-E 2 (-V,q)+*An/2). 
C\ + C 2 

where q = —ne + q is the charge of the island, the total r la n +i,n denotes transition from the left reservoir to the 

capacitance is \/C = 1/C\ + I/C2 and K% = C/Ci. The island, so that the number of electrons on the island is 

number of excess electrons on the island is n. The offset changed from n to n + 1, etc. and 

charge go is controlled by the gate voltage coupled to the p 

F{E) = 



island by the gate capacitance, q Q = C g V g . 1 exp ( E / k BT) 

At a finite current through the island, the energy di- where ksT is the thermal energy. We also define the total 

agram in Fig. ^ should be considered. Here we show the forward rates r l = r 1 ^ + F^ , where i — 1,2, and analo- 

equilibrium chemical potentials in the left and the right gous backward rates. The combined rate for an increase of 

reservoir and the non-equilibrium chemical potentials for the number of excess electrons on the island is r n+ i^ n — 

the spin-up and the spin-down electrons in the island. It is n +i,n + ^n+i.n an d analogous r n n+ i = r 1 n n+ i + 

assumed that the energy relaxation in the island is much r^ n , n+ i. The master equation which determines the prob- 
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ability p n of having n excess electrons on the island is of states. The spin-flip relaxation rate is 



dp n _ , r . p , 

—T- — —Pn [1 n-l,n + i n+l,n) 

dt 

+Pn+l^n,n+l + Pn- 1 -^n,n-l ■ 



(4) 



In the stationary state dp n /dt = and the detailed bal- 
ance symmetry gives r n+ i^ n p n = r ntn+ ip n+ i. The current 
through the first junction is 1\ = (/J + J^) , where the cur- 
rent of electrons with spin a is 



_ T-tler 
n+l,n ^ n — l.n 



In the second junction the current is I2 = (l\ + £ 



2)- 



T a — a n ( r 2 ^ r 2<T \ 

1l 2 — / y \ n—l,n 1 n+l,n j ■ 

n 

The transport of spins into the island is determined by 

out 



dt J 



tr 



dt 



{l\-ll-li+I l 2 )le. 



dt 



(5) 



In the stationary situation we can use the current conser- 
vation l\ + il = l\ + l\ and find 

The spin balance is 



ds / ds\ ( 

dt \dtJ tT \dt) rel ' 



(7) 



where (ds/dt) re j is the spin-flip relaxation rate. In equi- 
librium there are sq spins on the island (for a normal 
metal island sq — 0). The chemical potential is /io- The 
non-equilibrium spin-dependent chemical potentials are 
/i| = fio + 8/1 + A/i/2 and /ij = fio +S/i — A/i/2. The total 
number of spins on the island is s = so + [(S(i+Aii/2)5^~ 1 — 
(S/i — Afj,/2)57 ], where S^ 1 is the spin-dependent density 



ds\ 



A\i 



(8) 



.*7rel T sf T si S 
In the limit of fast spin-flip relaxation (r g £ — > 00), the 
number of spins on the island is s = sq + Sf^S^ 1 — S^ 1 ). 
Eq. (Q) determines the non-equilibrium chemical poten- 
tial shift A/i. Here T g £ is the spin-flip relaxation time in 
the island and <5 _1 = (S^ 1 + 8^ 1 )/2 is the average den- 
sity of states for spin up and spin down electrons at the 
Fermi level in the island (the inverse single-particle en- 
ergy spacing). Eq. (|^) effectively includes many-body ef- 
fects. The excess spin is related to the energy difference be- 
tween spin-up and spin-down states by s — s = XsA^i/ [i 2 B , 
where \s is the spin susceptibility (for noninteracting elec- 
trons Xs — f-B^ 1 ). In the stationary state (dp n /dt = 0, 
ds/dt = 0), the spin balance (^) can be written as 



L => [ = G s 2An/e, 



where the spin-flip conductance is defined by 

p 2 

G s 



(9) 



2tJ 



(10) 



From the equations for the tunneling rates (|^) and the 
spin balance (0), it is not directly obvious that there is 
only one unique solution for the non-equilibrium chemi- 
cal potential difference Afj, for a given set of parameters 
(temperature, gate voltage, source-drain voltage, capaci- 
tances and conductances). If Eq. ([?]) would be fulfilled for 
multiple solutions Afi for a given set of parameters, the 
current- voltage characteristics has a hysteretic behavior, 
with a solution Afi depending on the history. However, by 
extensive analytical and numerical studies with many dif- 
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ferent parameters to be presented below we always found Ferromagnet A) with magnetization configuration jf j or 
only one unique solution of Eq. (0). | J,f , which is the case for the device measured in Refs. ||, 

In this orthodox model the problem can be mapped on . Thus in this case our theory reduces to those in Refs. 
the equivalent circuit in Fig. [| by introducing a "spin-flip where the spins were assumed to be equilibrated, 

capacitance" C s = e 2 /2S, so that Second, we find that for a general system in the Coulomb 

blockade regime the current is zero, 1 = 0. and hence 

(es)/2 = C a (AfjL/e), Aft/ a = e 2 /(2G s ) =5. 

a vanishing non-equilibrium chemical potential difference 

This " charging energy" of the spin- flip capacitance is thus A „ . , ,. , , ^ , , , , , , 

° ° ° A/i = is a solution m the Coulomb blockade regime. 

simply the single-particle energy cost of a spin-flip, 6, or „, . . , . ,, . , , ,, 

° OJ ihis is physically reasonable since it means that there is 

more generally, the inverse of the magnetic susceptibility . , . . , , . , , , . , 

no spin accumulation on the island when tnerc is no cur- 



Vb/Xs 



rent flowing through the junctions. The Coulomb gap in 



Let us first discuss the situations in which An van- , . , . ,. n , , 

^ the current-voltage characteristics is not modified by the 

ishes identically in the stationary state. From the expres- . , , . 

spin accumulation. 

sions for the spin transport ([5]), the current conservation 
through the two junctions, I = I\ = I2, and the spin- 

dependent rates |), we see that when Af, = the spin- 3 Analytical results for a halfmetallic 
dependent currents are related to the total current by ferromagnet/normal metal/ halfmetallic 

t1 = t ferromagnet FSET 

1 1 + G u /G it ' 



l} = 



1 



I + Gii /Gil 



We will now discuss an idealized case for which an analytic 
where i — 1,2, expression for the current- voltage can be derived at zero 

j _ 2 t ( \ \ temperature. Let us assume that the leads are halfmetallic 

\1 + G11/G11 I + G2JG2]) 

ferromagnets and the island is a normal metal. The density 



From this relation we can make the following observa- 
tions. First, we see that in the symmetric case, G\-\ /Gij. = 
G21 /G21 , the spin-current is zero for A\x = and therefore 
A/i = is a solution of the spin balance equation (Jij) for 
any source-drain voltage V , gate voltage V g and temper- 
ature T, i.e. there is a solution with zero non- equilibrium 
chemical potential shift. This symmetry is expected for a 



of states at the Fermi energy vanishes for the minority 
spins in a halfmetallic ferromagnet. The capacitances are 
taken to be symmetric, C\ = C2 = C, and the gate voltage 
is tuned so that the offset charge is zero qo = 0. When the 
direction of the magnetizations in the ferromagnetic leads 
are antiparallel the tunnel conductances are 



structure of the type (Ferromagnet A- Ferromagnet B - G^ P = Gi, G^ P = 0, G^f = 0, G^\ P = G 2 . (11) 



8 Arne Brataas et al.: Non-equilibrium spin accumulation in ferromagnetic single-electron transistors 

The conductances G AP and G AP are zero, since there are where we have defined the conductance for the two junc- 

no spin-down states at the Fermi surface in the left reser- tions in series, Gyp — G\G 2 I iGx+G?). The non-equilibrium 

voir and there are no spin-up states at the Fermi surface chemical potential difference is determined by Eq. (||) which 

in the right reservoir. The electric current is I = l\ and is simplified to I = I s /2 = G s A[i/e as 
the spin-current is . G\ 2 / T ,AP 

^=orfeK - E < 

(ds \ T T 

~d~t J ^ ~~ ~ The current in the system can then be found to be 

a G G f E \ 

The spin-current is directly proportional to the current for / = s 1,2 j y m 1 j (13) 

G s + Gi,2 V e / 

any temperature T and source-drain voltage V. In the ab- 

at specific voltages determined by 

sence of spin-flip relaxation, r g £ — > 00, the current through 

eV AP =E (2rr Gs + Gl ' 2 I 1 
the system vanishes, because an electron cannot propagate m c \ G s 

from the left reservoir to the right reservoir without spin- We see from the current (|l|) that the Coulomb blockade 

flip. The master equation can be solved exactly at zero threshold is unchanged as expected. On entering the is- 

temperature for specific voltages V AP @: land from one of the leads the electrons must flip their 



a t> \ f 1 \ spins in order to be able to tunnel through the other lead. 

The total conductance of the system is thus given by the 

where E c = e 2 /2C is the Coulomb charging energy and m tW resistanceg in serieSj 1/G AP = 1/G> + 1/G± + ^ 

is an integral number. The non- vanishing; rates are u c ^ , . 1 • P \ 

to to as can be seen horn the equivalent circuit m iig. |3|. Mo 

Tut 1 n ri / \rM \ ; in \ instabilities (i.e. multiple solutions of the non-equilibrium 

r '„+!,« = -oGiE c {m-n)0{m-n), (12) M 



^-rr 1 chemical potential difference Afj, for a given set of param- 

F n,n+i = — ^GiE c (m - n)Q{n - m), 

< ^ eters) exist. 

r 2i n+ x,„ = — ^G 2 E c {m + n + l)<9(-m - n - 1), 

^ For a parallel orientation of the magnetizations in the 

r^n.n+i = — G 2 E c (m + n+ l)0(m + n + 1), 

e two leads, where A[i = 0, the current-voltage characteris- 

where 0{x) is the Heaviside function. The probability of tics is |p2f 

finding n electrons on the island is (n < m) I P = G\ :2 (V P ~)> (14) 

, \ ( G\\ n (ml) 2 . . where eV P = E r (2m + 1). Note that the specific volt- 

Pn(m) = — t [-7^7 — ,, m PoM 

\G 2 J (m — \n\)l(m + \n\)\ 

ages for which we have determined the current differ for 

and po (m) can be found from the normalization condition . . . 

the parallel and the antiparallel situation. The junction 

J2 n Pn( m ) = 1 H- The current is . . . 

' — p magnetoresistance is determined by 

7 = Gl - 2 ( y ™ P - ' JMR = (I P - I AP )/I P (15) 
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at the same voltages, which we cannot calculate exactly, current and the spin-current are calculated in units of the 

At high source-drain voltages at which the Coulomb charg- typical current Ge/2C. 

ing effects can be disregarded we see from Eq. (113) and T , c , . , , . T 

00 1 ■ — -J Let us first consider the DC transport properties. In a 

Eq. (114) and Fig. that the magnetoresistance for the , , ,. r , , , , , , , c 

1 v|_j/ o □ t> calculation of the current through the system for a given 

sys em is source-drain voltage V, we must first determine the non- 

G\ 2 

JMR = — '— — . equilibrium chemical potential difference A/i by the spin 

& s + &1,2 

In the limit of no spin- flip relaxation, r g f — > 00 (G s — > 0), 
the junction magnetoresistance is 100% since the current 
vanishes in the antiparallel configuration and in the limit 
of perfect spin-flip relaxation, r s £ — > (G s — > 00), the 
junction magnetoresistance vanishes. 



balance on the island (g). The solution of this equation is 
given by the intersection of the spin-current into or out 
of the island I s (Afj.) and the straight line G s 2Afi/e which 
describes the spin-flip relaxation in the island, where the 
spin-flip conductance G s was defined in Eq. (Eoh . 



4 General configurations: Numerical results 



As we have pointed out in the previous chapter, the 
non-equilibrium spin does not modify the Coulomb gap. 
However it does affect the current when the source-drain 



and discussions voltage is larger than the Coulomb gap. We show in Fig. 

|] the current as a function of the source-drain voltage V 

In the general situation with arbitrary junction conduc- f or a system consisting of G^/G = 0.3, G u /G = 0.1, 

tances and capacitances and at finite temperatures, the G21/G = 3, G2.1/G = 6, C\ = C2 = C, qo = and the 

spin-current and the electric current in the system have thermal energy k B T = 0.02E C . The upper curve shows 

to be calculated numerically. In experiments the tunnel the current when spin relaxation is fast G s /G = 1000, 

conductances of the two junctions depend strongly on the an d the lower curve shows the current in the case of slow 

thickness of the oxide tunnel barriers. We therefore present spin relaxation, G s /G = 0. The Coulomb blockade thresh- 

numerical results for a variety of possible realizations of ld for low source-drain voltages is clearly seen and the 

the junction parameters. steplikc structure for higher voltages called the Coulomb 

We choose symmetric capacitances C\ — C2 — C in staircase resulting from the discrete charging of the island 

our calculations. The important energy scale is then the M. The Coulomb blockade and the Coulomb staircase 

Coulomb energy E c = e 2 /2C and we scale the other rele- are smeared by the thermal fluctuations. Since the con- 

vant energies by this energy. The spin-dependent junction ductances of the tunnel junctions are spin-dependent, the 

conductances in the two junctions are described in units spin-flip relaxation on the island is important for the cur- 

of a typical junction conductance G, so that the electric rent. The current increases with increasing spin-flip con- 
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ductance G s (decreasing spin-flip relaxation time r s f). In with increasing spin-flip conductance (decreasing spin-flip 
the case of fast spin-flip relaxation [G s /G = 1000), the relaxation time ). As usual the peaks in the conductance 
spin-accumulation is vanishing small, and the steps occur as a function of the gate voltage increases with decreasing 
at eV = (2n + l)e 2 /2C, where n = 0, 1, 2, ... . For G s /G — temperature (not shown) Q. 
0, there is a spin- accumulation when eV > e 2 /2C, so that 

The junction magnetoresistance is the relative differ- 

the relative energy difference associated with the tunnel- 

ence in the resistance on switching the directions of the 

ing of one electron onto the island Ei(V,q) — aA/j,/2 is 

magnetizations in the leads from parallel to antiparallel 

spin-dependent. Hence, the spin-degenerate staircase at 

( |l5| ) . Let us consider the situation where the leads are fer- 

eV = (2n + l)e 2 /2C for n = 1,2,3... splits into two 

romagnetic and the island is non-magnetic, i.e. a F/N/F 

peaks, where the splitting is proportional to the spin- 
junction. The spin-accumulation causes a non-zero mag- 
accumulation. This can be seen in the lower curve in Fig. 

netoresistance. In the parallel configuration, the conduc- 
ts, where the second staircase (2n +1 = 3) occur at a 

tancesare G la = G 1 (l + aP)/2 and G 2(T = G 2 (l + crP)/2, 

voltage lower than 3e 2 /2C [^3| . 

where P is the polarization of the ferromagnet. We show 
Another typical experiment on the single-electron tran- in Fig. || the calculated junction magnetoresistance for 
sistor is to measure the differential conductance in the lin- G\ — G , G2 — 3G, V g — 0, fcsT = 0.05E C , C\ = C2 = C 
ear source-drain voltage regime as a function of the gate and a polarization P = 0.3 in the limit of slow spin- 
voltage V g (qo = VgCg). We show in Fig. || the influence flip relaxation G s = (upper curve), intermediate relax- 
of the spin- flip relaxation in the island on the differential ation G s /G = 1 (mid curve) and fast spin- flip relaxation 
conductance. The system is described by G-^/G = 3.0, G s — AG (lower curve). At low source-drain bias, we ob- 
G11/G = 0.2, G21/G = 0.1, G21/G = 2.0, Gi = G2 = G serve the magnetoresistance oscillations already reported 
and the thermal energy is fcsT = 0.1_E C . The curves show in |Tl],[l2|]. The amplitude of the oscillations decreases with 
the differential conductance as a function of the gate volt- increasing source-drain voltage, thus decreasing impor- 
age. The upper curve is calculated for fast spin relaxation tance of the Coulomb charging The period of the 

G s /G = 500, the mid curve is for intermediate spin re- oscillations for our system is close to 2E C . There is only 
laxation G s /G = 0.5 and the lower curve is for slow spin a small distortion of the shape of the magnetoresistance 
relaxation G s /G = 0.005. The typical oscillatory depen- oscillations with increasing spin-flip relaxation rate in the 
dence of the differential conductance with respect to the island. The oscillations in the TMR as a function of the 
gate voltage is seen. In the same way as the current at high source-drain voltage can be understood as a consequence 
voltages increases with increasing spin-flip relaxation in of the spin-accumulation in the antiparallel configuration. 
Fig. W, we see that the differential conductance increases The spin-accumulation increases with increasing current 
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through the system. We have seen in Fig. || that the cur- lem by numerically integrating Eq. (|J) and Eq. (^|). In 
rent has a steplike behavior with a period close to 2E C the upper part of Fig. (?] we show the calculated time- 
due to the discrete charging of the island. Hence the spin- dependent chemical potential difference when the source- 
accumulation also shows a steplike behavior (not shown) drain voltage is suddenly reduced at t — 0. The junctions 
as a function of the source-drain voltage. The magnetore- are described by G s /G — 0.5, G\^/G = 0.05, GiJG = 1, 
sistance for the F/N/F FSET increases with increasing G^jG = 2, G21/G = 0.01, G\ = C2 — C and the thermal 
spin-accumulation and hence shows oscillations with a pe- energy is ksT = 0.05-E c . We consider the high voltage case 
riod 2Ec- The magnetoresistance is noticeable even when V = 8E C , where the associated stationary electric current 
the spin-flip conductance is of the same order as the tunnel is I = 2.3Ge/2G and lower the source-drain voltage to 
conductances in agreement with Eq. (|l|). Disregarding the V — 2E C , where the associated stationary electric current 
Coulomb charging energy, the junction magnetoresistance is positive, I = 0.40Ge/2G. We have used i g f = 20RC, 
is which appears to be a reasonable estimate for junctions 

1 _ <v2 

JMR = P 2 - , (16) with a Coulomb charging energy of lOmeV and a junction 

1 — P 2 7 2 + a 2 ' 

conductance of R/Rk = 10 giving a charge relaxation 

where 7 = (Gi — G2)/ (Gi + G2) is a measure of the asym- 

time of RC = 0.5 • 10 _12 s. The chemical potential differ- 

metry of the junction conductances and a 2 = 4G S / (Gi + 

ence is seen to decay on a time scale much larger than 

G2) determines the reduction of the magnetoresistance 

the charge relaxation time. Finally we show in the lower 

due to the spin-flip relaxation. At a high source-drain 

part of Fig. |?] the current through the first junction (full 

bias, the numerical results agree well with Eq. (|l6|), giv- 

line) and the second junction (dotted line). It is clearly 

ing JMR= 6.9% for G s /G = 0, JMR= 3.4% for G s /G = 1 

seen that the relaxation of the current is slow on the time 

and JMR= 1.4% for G s /G = 4. 

scale RC. For time scales up to the order of RC, the cur- 

So far we have shown how the spin accumulation influ- 

rents through the first and the second junction differ due 

ences the DC transport properties. Spin accumulation also 

to the depopulation of the charge island. The Coulomb 

causes novel features in the AC and transient response of 

charging effect shows up as the almost constant current 

the system. Let us first fix the source-drain voltage at a 

on intermediate time-scales. 

high bias until the system is stationary and then lower 

the source-drain voltage. The transient current response In order to understand the dynamics it is useful to 
in this situation can be reversed on time scales of the inspect the device without the Coulomb charging effects, 
order of the spin-flip relaxation time, which is an unam- i.e. the capacitances Gi and G2 in the equivalent electric 
bigous signature of a non-equilibrium spin. In order to circuit in Fig. ^. We set the voltage on the left lead to zero 
explicitly show this we solved the time-dependent prob- and apply a time dependent potential V(t) to the right 
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lead. The complex impedance Z s ^ n (uj) — V(u>)/I(u>) is above (as used in Fig. ^). In this case the non-equilibrium 
1 G1G2 G11G2X — G11G21 A(j,(u>) ^ ^ spin accumulation decays slower since the spins must relax 



■^spin( w ) Gi + ^2 (^1 + eV(uj) through the spin- flip conductance G s on the island and the 

wnere transport through the junctions is suppressed. The spin- 

Afi(oj) _ 1 G11G2X — G*i|G , 2| . . accumulation time is then roughly equal to the spin- flip 

eV(to) ~ 1 + iwT spin {G a + G')(Gi +G 2 y 1 ' r 

relaxation time pCfl . The transient response to switching 

Here the spin accumulation time is 

on the source-drain voltage is similar. If the initial and 

Q 

t„ : = — , (19) final voltage are above the Coulomb charging energy, the 

" G s + G' 

where 1/G' = l/(Gi T +G 2T ) + l/(G u + G 2 |). From the re- 
lations @ and (£8p we see why switching-off the source- 
drain voltage (Vf = 0) reverses the transient current as 
found in the lower panel in Fig. Without the Coulomb 
blockade this transient decays on the time scale T g p m . In 
the limit that the junction conductances are much smaller 
than the spin conductance, the spin accumulation time 
( |l9| ) reduces to the spin-flip relaxation time, T g p m ~ r g f. 
In the opposite limit where the junction conductances are 
much larger than the spin conductance, the spin accumu- 

= /(^o)-JM f2Q) 

lation time is T S pj n ~ C S R. The spin-flip capacitance is ^ ' I(u> = 0) 
much larger than the charge-capacitance G in the regime as a function of frequency for an applied voltage which 
where the orthodox theory is valid (6 <C Eq) and thus the fluctuates around the Coulomb blockade threshold volt- 
spin accumulation time is much larger than the charge- age, V(t) — £? c (l+0.25 coswt). The temperature is kgT = 
relaxation time. The spin accumulation time in Fig. [7] 0.05E C , qo = 9 and C\ = Ci. The system is in the antipar- 
agree well with the value £ S p m = 0.43t g r as can be found allel configuration so that there is a spin-accumulation 
from the equivalent circuit neglecting the Coulomb block- on the island, Gij/G = 1.5, GiJG = 0.5, G^/G = 
ade. The spin accumulation time deviates from Eq. ( p^ ) if 0.5,G2x/G = 1.5,G S /G = 0.5 and the spin-flip relaxation 
the initial or final voltage is less than the Coulomb charg- time is r g r = 20RC. It is seen that the DC current varies 
ing energy pf]| . We show in Fig. ^ the average current most strongly when wr g f ~ 1, which is a consequence 
through the first and the second junction when the fi- of the spin-dynamics in the system with the character- 
nal voltage is zero with the same system parameters as istic time-scale r g £. A corresponding calculation with the 



response time is determined by Eq. (19) otherwise the re- 
sponse time roughly equals the spin-flip relaxation time. 

The long time response of the system due to the spin 
dynamics can also be observed in other AC transport ex- 
periments. A fast single-electron transistor has recently 
been realized p4|. We therefore study the influence of an 
AC source-drain voltage V(t) = Vq + V\ cosojt on the DC 
current through the system, I(uo) = (j/(2tt) J^° +2 ^^ dtl(t), 
where to is an arbitrary time constant. We show in Fig. |^ 
the relative change in the DC current 
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same parameters in the parallel configuration shows no members for their hospitality at the NTT Basic Research Lab- 
frequency dependence of R(u) around CJT g r ~ 1. oratories and Keiji Ono for a discussion. 
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tance increases with increasing spin-flip relaxation rate. 
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are time-dependent. 
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Fig. 5. The conductance dl/dV as a function of the gate- 
voltage V g for different spin-relaxation conductances: G s /G = 
500 (upper line), G s /G = 0.5 (mid line) and G s /G = 

0. 005 (lower line). The system parameters are Gif/G = 3.0, 
G n /G = 0.2, G 2T /G = 0.1 , G 2i = 2.0, T = 0.1E C , go = 0, 
Gi = C 2 = 1.0. 

Fig. 6. The magnetoresistance as a function of the source- 
drain voltage in a F/N/F double junction system with P=0.3, 
Gif + Gi| = G, G 2 | + G 2 | = 3G, V g = 0, Ci = G 2 = 1 
and ksT — 1.0. The upper curve is for slow spin-relaxation, 
Gs/G — 0, the mid curve for intermediate spin-flip relax- 
ation Gs/G — 1 and the lower curve for faster spin- relaxation, 
Gs/G = 4. 

Fig. 7. The current as a function of time. The source-drain 
voltage is switched from V = 8E C to V = 2E C at t = 0. The 
system parameters are G s /G = 0.5, Gif/G = 0.05, G14./G = 

1, G 2T /G = 2, G 2J ./G = 0.01, k B T = 0.05E C , go = and 
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Fig. 8. The average current through the first and the sec- 
ond junction as a function of time. The source-drain voltage is 
switched from V = SE C to V = at t = 0. The system param- 
eters are G s /G = 0.5, Gi T /G = 0.05, G u /G = 1, G 2T /G = 2, 
G 2J ./G = 0.01, k B T = 0.05£ c , q = and Ci = G 2 

Fig. 9. The time-averaged (DC) current as a function of the 
frequency of the applied source-drain voltage, V(t) — Ec(l + 
0.25 cosuit). The reservoirs are in the antiparallel configuration, 
Gi T /G = 1.5, G u /G = 0.5, G 2T /G = 0.5,G 2J ./G = 1.5. The 
other parameters are fcsT = 0.05i5 c , qo = 9 and C\ = C2 
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